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Abstract 

The a decay half-lives of recently synthesized superheavy nuclei (SHN) are calculated 
by applying a new approach which estimates them with the help of their neighbors 
based on some simple formulas. The estimated half-life values are in very good agree- 
ment with the experimental ones, indicating the reliability of the experimental ob- 
servations and measurements to a large extent as well as the predictive power of our 
approach. The second part of this work is to test the applicability of the Wentzel- 
Kramers-Brillouin (WKB) approximation for the quantum mechanical tunneling prob- 
ability. We calculated the accurate barrier penetrability for alpha decay along with 
proton and cluster radioactivity by numerically solving Schrodinger equation. The 
calculated results are compared with those of the WKB method to find that WKB 
approximation works well for the three physically analogical decay modes. 

PACS: 27.90. +b, 21.10.Tg, 23.60.+e, 23.50,+z, 23.70.+j 
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1 Introduction 



Over the past decades, the syntheses of superheavy elements and their lifetime measurement 
have been explored with a variety of methods. The heavy elements with Z = 107 — 112 
have been successfully synthesized at GSI [lj. Elements along with Z = 113 — 116, 118 have 
been produced at JINR-FLNR, Dubna [2]. Last year, two isotopes of a new element with 
atomic number Z = 117 were synthesized in the fusion reactions between 48 Ca projectiles 
and radioactive 249 Bk target nuclei whose a chains terminated by spontaneous fission was 
observed in Dubna [3J, which fills the gap between the elements 116 and 118. The element 
114 was independently confirmed recently by the LBNL in the USA [4] and GSI [5]. A 
superheavy element isotope 285 114 was observed in LBNL last year [6], and an isotope of 
Z = 113 has been identified at RIKEN, Japan [7j. Thus up to now superheavy elements 
with Z = 104 — 118 have been synthesized in experiment and consequently they offer the 
possibility to study the heaviest known nuclei with greater detail. However, their is no 
consensus among theorists with regard to what should be the next doubly magic nucleus 
beyond 208 Pb. Nearly all of modern calculations predict the existence of a closed neutron 
shell at N = 184. However, they differ in predicting the atomic number of the closed 
proton shell. For instance, the macroscopic-microscopic model predicts the shell gap at 
Z = 114 [8TfT0]. The microscopic Skyrme-Hartree-Fock models give Z = 124,126 [TTHT3] 
and the relativistic mean-field calculations suggest Z = 120 [T4HT6] . The magic numbers 
Z = 132 and iV = 194 were predicted from the discontinuity of the volume integral at shell 
closures [17J. A tremendous progress in experiments and the development of the radioactive 
ion beam facilities have made it possible to reach the island of superheavy elements. 

The heaviest SHN decay primarily by the emission of a-particle terminated by spon- 
taneous fission. Therefore, in recent experiments, a decay has been indispensable for the 
identification of new nuclides. Because the experimentalists have to evaluate the values of 
the a decay half-lives, during the experimental design, it is quite important and necessary to 
investigate the a decay of SHN theoretically. Although a decay is very useful for the study 
of the nuclei, a quantitative description of them with a satisfying accuracy is difficult. The 
a decay was firstly interpreted as a consequence of quantum penetration of a-particle by 
Gamow in 1928. At present, many theoretical approaches have been being used to describe 
the a decay, such as the cluster model [T8T[20] . the density-dependent M3Y (DDM3Y) effec- 
tive interaction [2"T||2"2"]. the generalized liquid drop model (GLDM) [2"3Tf2"T] . the Coulomb and 



proximity potential model [28J, the superasymmetric fission model [29||30]. the UMADAC 
method [31 J, the coupled channel approach [32|l33] and the universal curves for a and cluster 
radioactivities in a fission theory [33] . Some physically plausible formulas also were employed 
to calculate the a decay half-lives directly [29 |l35H39] . Interestingly, the superasymmetric 
fission theory for a and cluster decay has been extended by some authors to study metallic 
cluster physics [40J[4T] , which is an example of using the nuclear methods in nanophysics. 

The half-life is extremely sensitive to the a decay Q value and an uncertainty of 1 MeV 
in Q value corresponds to an uncertainty of a-decay half-life ranging from 10 3 to 10 5 times 
in the heavy element region [32] ■ In this work, with the experimental Q values, we carry 
out the half-life calculations for the recently synthesized SHN by employing a relationship 
between the a decay half-lives of neighboring SHN that are established based on some 
simple semi-empirical formulas. Differently from our approach, theoretical estimates for the 
lifetimes by calculating the quantum mechanical tunneling probability in a WKB framework 
is widely performed for the a decay along with other physically analogical decay processes. 
It is pointed out that the WKB approximation works well at energies well below the barrier 
height [33]. As a matter of fact, the accuracy of the WKB approximation also depends on 
the shape of the potential barrier as well as the decay energy. In this work, we obtain the 
penetrability with a different method and show the applicability of the WKB approximation 
in a, proton and cluster radioactivity. 

This paper is organized as follows: In Section 2, a brief discussion of the method and 
the calculated results along with the corresponding discussions for the half-lives of SHN are 
presented. The applicability of the WKB approximation for a decay, proton and cluster 
emission are discussed in Section 3. Finally, a brief summary is provided in Section 4. 



2 a decay half-lives of superheavy nuclei within a new 
approach 

We start from Royer's [44J and Viola-Seaborg semi-empirical (VSS) formulas j3HJ[36]. The 
Royer's formula is given by 

cZ 

log 10 T(s) = a + L4 1 / 6 v / Z+^=, (2.1) 



where Q is in MeV and the parameter set varies for four types: a = —25.31, b = —1.1629, 
c = 1.5864 for even(Z)-even(iV), a = -26.65, b = -1.0859, c = 1.5848 for even-odd, 
a = -25.68, b = -1.1423, c = 1.592 for odd-even and a = -29.48, b = -1.113, c = 1.6971 
for odd-odd nuclei [44J. In our previous work, this formula has been extended by taking 
into account the centrifugal barrier to describe unfavored a decay. For odd-mass nuclei, it is 
possible that some decays involve non-zero I values. However, as no experimental evidence 
is available for the spin-parity of the levels involved in the decay, we have not included the 
centrifugal barrier in the present calculations. The VSS formula is given by 

log 10 T(s) = {aZ + b)^= + cZ + d + h log . (2.2) 

v Q 



Instead of using the original set of constants by Viola and Seaborg, recent values a = 1.64062, 
b = —8.54399, c = —0.19430, d = —33.9054 being valid for the nuclei of four types are 
used [37]. h\ og accounts for the hindrances associated with odd proton and odd neutron 
numbers but does not take an effect in our calculations. Once the half- life of a nucleus A °Z 
(reference nucleus) is known, the half-life of an other nucleus A Z (target nucleus) with the 
same type can be derived. The difference of the logarithms of half-life is written with Eq. 
M) as 

S = log 10 T — log 10 T 

. \ / 7, Zr. \ 

(2.3) 



= b^VZ-A^^j+cf 
and with Eq. ( 12.2ft as 



VQ VQ~o 



S = (aZ + 6) J= - (aZ + b)^= + c(Z - Z ). (2.4) 
vv vvo 



Therefore, the half-life of the nuclei Z can be obtained from T = 10 S T with the help 
of its neighboring nucleus A °Z . The two formulas can validate each other to obtain more 
compelling results. 

Now we focus on two simple cases. One is that the two nuclei are in an isotope chain for 
which Eqs. (12. 3p and ( 12. 4p can be further simplified, and the other is that the two nuclei 
belong to an a decay chain. We estimated the a decay half-lives of recently synthesized 
SHN A Z with the help of the reference nuclei A ~ 2 Z and A+2 Z by employing Eqs. (12. 3 j) 
and (12.41) without taking into account the uncertainty of the experimental Q values. The 
results are presented in Table 1 compared with experimental data. The results obtained 
with the DDM3Y effective interaction and the GLDM are also shown for comparison. The 



third column marks the experimental a decay half-lives, and the columns 4 and 5 are the 
estimated ones from Eq. (12 .3p based on the Royer's formula and from Eq. (12. 4p based on 
the VSS formula, respectively. The first half is obtained with the reference nuclei A ~ 2 Z and 
the second half are obtained with the reference nuclei +2 Z. On an average, the DDM3Y 
results are slightly larger than the experimental data while the GLDM values are lower than 
the measured ones. In fact, the DDM3Y effective interaction and GLDM are very successful 
because of the appropriate treatment on the microscopic level in the DDM3Y calculation 
and the quasimolecular shape in the GLDM in consideration of the difficulty in accurate 
half-life calculation. Our calculated values are in very good agreement with the experimen- 
tal measurements which indicates our method is a very effective approach to investigate the 
half-lives of a decay when the experimental Q values are given though it is very simple in 
theoretical framework compared to the DDM3Y and GLDM. The two approaches based on 
the Royer's formula and VSS formula give nearly the same results implying the reliability 
of our method to a certain extent. The half-life of 282 113 is underestimated by one order 
of magnitude with the present method and the GLDM, which are possibly due to nonzero 
angular momentum transfer or some nuclear structure effects such as the dramatic deforma- 
tion changes as suggested in Ref. |48j and the influence of a possible neutron shell gap at 
N = 166 on its daughter nucleus. In Ref. [H], it is suggested that N = 166 is a neutron shell 
gap in certain region within relativistic mean field models. This nuclide warrants further 
experimental measurements with higher statistics. 

Apart from calculating the decay half-life, the present method is also a useful approach 
to validate the experimental measurements. The recently observed SHN still await confir- 
mation by other laboratories, which is not easy because the new SHN form an isolated island 
that is not linked through a decay chain with known nuclei. Therefore, the theoretical con- 
firmations become important and necessary. Since the half-lives of the reference nuclei are 
taken from the experimental values, that the estimated results with Eqs. (12. 3p and (12 .4p are 
in excellent agreement with the experimental values suggests that the experimental half-lives 
are themselves consistent with each other which confirms the reliability of the experimental 
observations and measurements to a great extent. 

In order to further confirm the above conclusions drawn from Table 1 and the predictive 
power of Eqs. (12 .3p and (12. 4p . a relationship between the a decay half-lives of the nuclei 
belonging to an a decay chain is investigated here and the estimated half-lives of SHN A Z 



are listed in Table 2. The columns 4 and 5 are the estimated half-lives from Eqs. ( 12. 3 j) 
and (12. 4p with the help of their daughter nuclei A ~ 4 (Z — 2) while the columns 6 and 7 
are the estimated ones with the help of their parent nuclei A+4 (Z + 2) respectively. The 
well agreement between the evaluated values with the experimental ones further indicates 
the predictive power of Eqs. (12.31) and (12. 4p along with the reliability of the experimental 
measurements. Of course, the uncertainties in the measured values are large because of the 
experimental difficulties and poor statistics. And our approach underestimates the a decay 
half-life of 279 111 as the DDM3Y and GLDM. The reason for that is just the same as that 
for 282 113 mentioned above. Although the main shell effect has been included in the Q value 
to evaluate the half-life, the preformation probability is also affected obviously by the shell 
effect [25], which will lead to a large deviation of Eqs. ( 12. ip and (12.21) together with Eqs. 
(12.31) and (12. 4p for the nuclei around the magic numbers. However, as shown in Table 2, the 
half-lives of element 114 and the isotones iV = 172 are very well reproduced by applying 
Eqs. ( 12.31) and ( 12. 4p . In other words, it does not exhibit any evidence to show that Z = 114 
and iV = 172 are shell closures in this region. 

Finally, we predict the a decay half-lives of the nuclei belonging to an a decay chain 
starting from 293 117 without experimental values by employing Eqs. H 2 . 3 1) and ( 12 .4p . They are 
listed in Table 3 which may be useful for future experimental measurements. The theoretical 
half-lives from Eq. (12. 3p and Eq. (12. 4p agree very well with each other which additionally 
confirms again the validity of our approach. The deviation in the predicted values for 290 115 
may be large since the experimental Q value has a large uncertainty of 0.41 MeV. The Q 
value for this nuclide needs to be measured with a higher accuracy. For other still unknown 
SHN, one can make predictions for the a decay half-lives with the decay energies calculated 
from the atomic mass evaluation of Audi et al. [SD] as a substitute since the agreement with 
the experimental data on the mass of the known heaviest elements is very satisfactory, or 
from a formula for the a-decay energy |51j . 

3 Applicability of the WKB approximation 

We turn now to the applicability of the WKB approximation for a decay. The interaction 
potential V(r) is the sum of the nuclear potential Vjv(r), Coulomb potential Vc(r) and the 
centrifugal barrier. We approximate the nuclear potential by Vjv(r) = — Ai?7o/[l+exp( r ~ Q Ro )], 



with #0 = 1.27A 1 / 3 fm, a = 0.67 fm, U = [53 - 33(iV - Z)/A] MeV, A x the mass number of 
the emitted particle, N, Z, A the neutron, proton, and mass numbers of the parent nucleus 
respectively. Here the single particle potential is taken from Ref. (52]. The Coulomb potential 
is given by the point-like plus uniformly charged sphere method with a parent nucleus radius 
R = 1.28A 1/3 - 0.76 + 0.8A~ 1/3 fm [53]. As an example, the barrier for 212 Po-^ 208 Pb+a is 
shown in Figure 1(a). Here only the barrier is considered and we divide the barrier into a 
sequence of square barriers, as shown in Figure 1(b). In principle, the barrier ranging from 
r\ to infinity should be taken into account for the calculation yet it is unpractical. Therefore, 
we cut off the barrier at a sufficiently large distance of r 2 = 1000 fm and the potential barrier 
is divided into n = 60000 parts with a step of h = (r 2 — r-y)/n without loss of accuracy. The 
wave function u{r) (^("F") = Yi m (9, <p)u(r)/r) of the emitted particle with Q value in these 
n regions can be written as 

u\ = A ltl exp(ik 1 r 1 ) + A 2 ,iexp(-ik 1 r 1 ), 
u 2 = A 1)2 exp(ik 2 r 2 ) + A 2)2 exp(— ik 2 r 2 ), 

u n _! = A lin _iexp(iA; n _ir n _ 1 ) + A 2>n _ 1 exp(-zA; n _ 1 r n _i), 
u n = A 1)n exp(ik n r n ) + A 2tn exp(-ik n r n ), (3.5) 



with kj = v/2/i(Q - Vj)/h 2 , r 3 = n + (j - l)h and Vj = [V(rj) + V(r j+1 )]/2. The wave 
function outside of the barrier is given by 

u = A lfi exp(ik r ) + A 2fi exp(-ik r ), 
u n+ i = A l!n+1 exp(ik n+1 r n+1 ), (3.6) 

with k = k n+ i = y/2/jQ/H 2 where A\$ and Ai >n+1 are the amplitude of incident wave and 
transmitted wave, respectively. By using the connection condition of wave function, one can 
deduce the transmission amplitude and reflection amplitude for the nth square barrier 

1/- k n+ \ 



A x ,n = -yl + -j—)exp[i(k n+1 -k n )r n }A lin+1 , (3.7) 
1 ( ^ k n+ \ 



A 2 ,n = ^ ( 1 ~ I eX P (*n+l + kn ) r ^ ^1."+!' ( 3 -^ 



and for the jth (j < n) square barrier 



A X j = - exp (-ikjTj) 



kj+i- 

(3.9) 



exp (ik j+1 rj) (1 + -j^-)A 1J+1 + exp (-ik j+1 rj) (1 - -t±)A ijj+1 



A 2 j = - exp {ikjTj) exp (ik j+1 rj) (1 J -^-)A 1J+1 + exp (-ik j+1 rj) (1 + -j^-)A 2J+1 



(3.10) 



The penetration probability is given by 



.4 
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l,n+l 



(3.11) 
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Normalization won't help-this is not a normalizable state. We choose Ai^ n+ i = 1, and 
then the v4 10 can be recured according to the above formulas and hence one can obtain 
the penetrability. As a matter of fact, our method is a kind of numerical method to solve 
one-dimensional Schrodinger equation for unbound state, in which the differential equation 
is translated to recursion formulas. 

Before we perform the calculation for a, proton and cluster radioactivity, we have checked 
this method with a soluble example 



for which the exact analytic transmission probability is known [5ljl55]. It is found that 
the result with this method completely coincides with the analytic one, which confirms the 
reliability of this method. Taking this method as a standard, one can test whether the WKB 
approximation works well or not. In the WKB approximation, the formula 



is employed to evaluate the penetrability, and one can estimate the relative deviation RD = 
(Pwkb — P)/P x 100% of this WKB method. As a semiclassical approximation, there exist 
two classical turning points r in and r out in WKB method. The penetration only performs 
between r in and r out and the effects of the regions I and III of potential barrier in Figure 1(a) 
are neglected. However, according to quantum mechanics, the particle can be also reflected 
back in the region III with some probability. Additionally, the WKB method also brings 
some errors when one evaluates the penetrability from r in to r out (the region II) and it cannot 
deal with these Q values being near or larger than the top of the barriers in principle while 
our fully quantum mechanical approach has no such a drawback. 

We select a decay events with 52 ^ Z ^ 118, and the experimental Q values are taken 
from Refs. |2][53]. For a decay, the RD values have been presented in Figure 2. The WKB 
approximation underestimates the penetration probability by about —40% ~ —30%. It 



V{x) = V cosh- 2 (x/a),V > 



(3.12) 




(3.13) 



is not possible to calculate the a decay half-life theoretically with a high accuracy within 
the framework of barrier penetration because the preformation factor is very difficult to be 
estimated microscopically and the a-daughter interaction has not yet been well determined. 
From this point of view, the WKB approximation thus works well in investigations of a decay 
especially for SHN since the experimental half-lives of SHN tend to have a large uncertainty. 
Indeed, because the deviation is nearly a constant in the whole mass region, this constant 
error can be compensated by other quantities such as a phenomenological assault frequency 
in actual calculations within the WKB framework. In an analogous way, we investigate 
the deviation of the WKB approximation for proton and cluster radioactivity. The study 
of the nuclei far away from the /3-stable line has attracted world wide attention from both 
the experimental and theoretical points of view. In the case of very proton-rich nuclei, it is 
expected to observe the proton emission experimentally [56]. Since around 1980, the cluster 
radioactivity was observed in experiments with daughter nuclei being almost closed-shell 
spherical nuclei around 208 Pb. The proton and cluster emission can be treated as simple 
quantum tunneling effects through a potential barrier just as the a decay [57JEE]- We 
select cluster emitters with emitted particles from 14 C to 34 Si and spherical proton emitters, 
for which the experimental Q values are taken from Refs. [501 ISHl ISH] - Figure 3(a) shows 
the relative deviation RD of the WKB method for the proton radioactivity of spherical 
proton emitters. As can be seen, the WKB method underestimates the penetrability by 
—40% ~ —20% and again the deviation does not fluctuate with a large amplitude as that 
for a decay. Figure 3(b) presents the relative deviation RD of the WKB method for cluster 
radioactivity. It indicates the WKB approximation works well for the cluster emission with 
a deviation by only —5% ~ 15%. The RD is found to be insensitive to the nuclear potential 
Vn{t) for these three decay modes which indicates the conclusion we draw here is universal. 



4 Summary 

To summarize, the a decay half-lives of newly synthesized SHN have been investigated 
in terms of the correlation between the half-lives of a decay. The results of the present 
calculations with this relationship based on the Royer's and VSS formulas are in excellent 
agreement with the experimental data which indicates the predictive power of our method. 
According to the present calculations, an important conclusion is that the experimental half- 



lives are themselves consistent with each other confirming the reliability of the experimental 
observations and measurements to a great extent. For the nuclei 282 113 and 279 111, the 
half-lives from the theoretical estimations are underestimated by one order of magnitude 
possibly due to nonzero angular momentum transfer or some nuclear structure effects. The 
half-lives of the synthesized SHN does not exhibit any evidence to show that Z = 114 
and N = 172 are shell closures in the considered region according to our analysis. The 
other task of the present work was to test the applicability of the WKB approximation for 
quantum mechanical tunneling probabilities. We calculated the barrier penetrability for a 
decay, proton and cluster emission accurately with the recursion formulas by dividing the 
potential barrier into a sequence of square barriers, which is a numerical method to solve one- 
dimensional Schrodinger equation for unbound state, and the results are compared with those 
of the WKB approximation. It is found that the WKB method produces relative deviations 
by about —40% ~ —30% for a decay of heavy and superheavy nuclei, —40% ~ —20% for 
proton emission and —5% ~ 15% for cluster radioactivity. Also, in consideration of the 
deviations being nearly constants in each decay mode, indeed, the WKB approximation 
works well for these decays. 
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Table 1: Calculated a-decay half-lives of recently synthesized SHN A Z within Eqs. (12.31) and 
(12.41) taking the nuclei in the isotope chains as references. The first eight nuclei and the rest 
ones are obtained with the reference nuclei A ~ 2 Z and A+2 Z, respectively. The experimental 
data [2] and other theoretical results are also listed for comparison. Some experimental Q 
values are obtained by using the measured a kinetic energies taking into account the electron 
shielding corrections. 



Nucleus Q expt - T cxpt - 



T |2~4l 



53^ ms 62.0l£$ ms 66.1 
18^ 6 ms 21. 2 +a5 
22+ - 26 s 

2-7_o.7 s 
0.48±g-f? s 
34^ 7 s 
3.6 ±tl s 
18±f ms 

7.ljt? ms 
32l}f ms 

o-48±S:Sg s 

73l^ 4 ms 
3.8iJ;? s 
4.2lTj ms 



22.0 

22.7^° 4 2 

u - lo -0.06 
1 7+0.58 
x - ' -0.33 b 

4.3^- 9 7 s 

ro c+16.6 

. U g y 

1 Q+3.4 
L - J -0.8 b 

14.4l* 6 2 9 
5.615 9' ms 

r 9 r>+61.7 

uz.z_ 19 q 
1+0.39 



yDDM3Y |22j J 1 



iGLDM 



206f^"ms 
39l 2 ° ms 



22.81+^ u f ms 
10.45±|;^i ms 



293 n6 
292 n6 
289 n5 
289 n4 
284 113 
285 n2 
280 m 
291 n6 
290 n6 

287 115 
287 n4 

282 113 

283 n2 

278 ln 



10.67(6) 

10.80(7) 

10.50(9) 

9.96(6) 

10.15(6) 

9.29(6) 

9.87(6) 

10.89(7) 

11.00(8) 

10.74(9) 

10.16(6) 

10.83(8) 

9.67(6) 

10.89(8) 



ms 



-5.1 
13+ - 61 s 

h0.54 „ 



ms 
ms 



0.31 



15.9 
-9.3 
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1.6" 

5.6i^ 3 
50.2 

2.9 1 2 >~> 

15.4^ 5 ° 
6.0±H ms 

e+65.5 
oo.o_ 202 

0.79l££ s 

6.31^ ms 

2.6±i 3 7 s 
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-1.9 



-1.8 
-1.2 b 
;+0.72 
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75+ 41 « 
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1.55" 



0.52"_ 
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ms 



1 q+0-9 
60.4 H 
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0.74l5;ig s 
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-26 
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13-4l 5 7 i 

n 7+35.8 
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o 4 7+1.99 
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u - iu -0.05 b 
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0-95l !32 s 
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Table 2: Calculated a-decay half-lives of recently synthesized SHN A Z within Eqs. ( 12.31) 
and (12.41) taking the nuclei in a-decay chains as references. The results in columns 4 and 5 
are obtained with the help of their daughter nuclei A ~ 4 (Z — 2) while the columns 6 and 7 
are the estimated ones with the help of their parent nuclei A+4 (Z + 2). The experimental 
data [2] and other theoretical results are also listed for comparison. Some experimental Q 
values are obtained using the measured a kinetic energies taking into account the electron 



shielding corrections. 
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Table 3: Predicted a-decay half-lives of the nuclei A Z in a-decay chain starting from 294 117 
within Eqs. (12.31) and (12 M . The experimental Q values are obtained with the measured 
a kinetic energies [3] taking account of the electron shielding corrections. The results in 
columns 3 and 4 are obtained with their isotopes A ~ 2 Z and the ones in columns 5 and 6 are 

obtained with the help of 294 117. 

nuclei Q Ex P(MeV) T^B T M yES 
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Figure 1: (a) Potential barrier of 212 Po— > 208 Pb+a. (b) It is divided into a sequence of square 
barriers. 
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Figure 2: Relative deviation of penetrability caused by the WKB approximation for a decay. 
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Figure 3: Relative deviation of penetrability caused by the WKB approximation for proton 
and cluster radioactivity. 



